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l> I ABSTRACT 

■ An analytical model is presented for the postcollapse equilibrium structure 

of virialized objects which condense out of a low-density cosmological back- 
] ground universe, either matter-dominated or flat with a cosmological constant. 

. This generalizes the model we derived previously for an Einstein-de Sitter 

(EdS) universe. The model is based upon the assumption that cosmological 



o 

?— j ' haloes form from the collapse and virialization of "top-hat" density perturba- 

^ | tions and are spherical, isotropic, and isothermal. This leads to the prediction 

of a unique, nonsingular, truncated isothermal sphere (TIS), a particular so- 
• i— ( . 

■ lution of the Lane-Emden equation (suitably modified when A^0). The size 

$H ' 

and virial temperature are unique functions of the mass and redshift of for- 
mation of the object for a given background universe. The central density is 
roughly proportional to the critical density of the universe at the epoch of 
collapse. This TIS model is in good agreement with observations of the inter- 
nal structure of dark matter-dominated haloes on scales ranging from dwarf 
galaxies to X-ray clusters. It also reproduces many of the average properties 
of haloes in simulations of the Cold Dark Matter (CDM) model to good accu- 
racy, suggesting that it is a useful analytical approximation for haloes which 
form from realistic initial conditions. Our TIS model matches the density 
profiles of haloes in CDM N-body simulations outside the innermost region, 
while avoiding the steep central cusp of the latter which is in apparent con- 
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flict with observations. The TIS model may also be relevant to nonstandard 
CDM models, like self-interacting dark matter, recently proposed to resolve 
this conflict. 

Key words: cosmology: theory - dark matter - galaxies: clusters: general - 
galaxies: formation - galaxies: haloes - galaxies: kinematics and dynamics 

1 INTRODUCTION 

Galaxies and clusters of galaxies formed when gravitational instability amplified density 
fluctuations in the expanding cosmological background universe. Regions denser than aver- 
age eventually stopped expanding and recollapsed to form virialized objects. The question 
of what equilibrium structure results when a density perturbation collapses out of the ex- 
panding background universe and virializes is central to the theory of galaxy formation. 
The nonlinear outcome of the growth of Gaussian-random-noise cosmological density fluc- 
tuations due to gravitational instability in a hierarchical clustering model like CDM is not 
amenable to direct analytical solution, however. Instead, numerical simulations are required. 
As a guide to understanding these simulations, as a check on their accuracy, and as a means 
of extrapolating from simulation results of limited dynamic range, analytical approxima- 
tions are nevertheless an essential tool. One such tool of great utility has been the solution 
of the spherical top-hat perturbation problem (e.g. Peebles 1980, Padmanabhan 1993). As 
used in the Press-Schechter ("PS") approximation (Press & Schechter 1974) and its various 
refinements, the top-hat model serves to predict well the number density of virialized haloes 
of different mass which form at different epochs in N-body simulations. An analytical model 
for the internal structure (e.g. mass profile, temperature, velocity dispersion, radius) of these 
virialized haloes would be a further tool of great value for the semi-analytical modelling of 
galaxy and cluster formation, therefore. Shapiro, Iliev & Raga (1998; Paper I) derived such 
a model for an EdS universe. Here we shall generalize the analysis of Paper I to the case of a 
low-density universe (Q < 1) which is either matter-dominated or flat with a cosmological 
constant (i.e. f2 = 1 — A ). We shall also describe how the TIS model can be generalized to 
other cosmologies with a nonclumping background component, such as quintessence. 

It is generally assumed that the collapse to infinite density predicted by the exact, non- 
linear solution of the spherical top-hat perturbation problem is interrupted by a violent 
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relaxation to virial equilibrium at a finite density, as a result of the growth of initially 
small-amplitude inhomogeneities in the density distribution. Earlier work adopted the crude 
approximation that the postcollapse object which emerges from this violent relaxation is ei- 
ther a uniform sphere or a singular isothermal sphere, with the same total mass and energy 
as the collapsing top-hat and with a radius and velocity dispersion (or, equivalently, gas 
temperature) fixed in accordance with the virial theorem. Our first motivation, therefore, 
was simply to improve upon this earlier treatment by finding a more realistic outcome for 
the top-hat problem. As a starting point, we adopted the assumption that the final equilib- 
rium is spherical, isotropic, and isothermal, a reasonable first approximation to the N-body 
and gasdynamic simulation results of the CDM model. The postcollapse analytical solution 
we derived in Paper I from this assumption quantitatively reproduces many of the average 
properties of the haloes found in those simulations to good accuracy, so we are encouraged 
to believe that our approximation is well justified for haloes which form from realistic initial 
conditions. Our model is in disagreement, however, with one detail of the N-body simula- 
tion results that, in their very centers, simulated haloes have cuspy profiles (e.g. Navarro, 
Frenk & White 1997; "NFW"). By contrast, our model predicts a density profile with a 
small, uniform-density core, in good agreement with the observed rotation curves of dark 
matter-dominated galaxies (Iliev and Shapiro 2001) and with cluster mass profiles inferred 
from gravitational lensing (Iliev and Shapiro 2000, Shapiro and Iliev 2000). This apparent 
discrepancy between the cuspy profiles of the N-body results and the observed dark matter- 
dominated haloes (e.g. Moore et al. 1999) has led recently to a vigorous reexamination of 
the cold, collisionless nature of CDM, itself, and the suggestion that a variation of the mi- 
crophysical properties of the dark matter might make it more "collisional" , enabling it to 
relax dynamically inside these haloes so as to eliminate the central cusp (e.g. Spergel & 
Steinhardt 2000; Burkert 2000; Dave et al. 2000, Firmani et al. 2000 and references therein). 
While the details of this suggestion are still uncertain, our model can also be applied in that 
case, to the extent that we are able to ignore the details of the relaxation process inside the 
halo and approximate the final equilibrium as isothermal. 

As described in Paper I, our model is as follows: An initial top-hat density perturba- 
tion collapses and virializes, which leads to a truncated nonsingular isothermal sphere in 
hydrostatic equilibrium (TIS), a solution of the isothermal Lane-Emden equation. Although 
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the mass and total energy of the top-hat are conserved thru collapse and virialization, and 
the postcollapse temperature is set by the virial theorem (including the effect of a finite 
boundary pressure), the solution is not uniquely determined by these requirements alone. 
In order to find a unique solution, some additional information is required. We adopted the 
anzatz that the solution selected by nature is the "minimum-energy solution" such that the 
boundary pressure is that for which the conserved top-hat energy is the minimum possible 
for an isothermal sphere of fixed mass within a finite truncation radius. This assumption 
fixes the ratio of the radius of the postcollapse sphere to the radius of the top-hat perturba- 
tion which created it, as measured at the latter's epoch of maximum expansion, uniquely. 
For comparison, we appealed to the details of the exact, self-similar, spherical, cosmological 
infall solution of Bertschinger (1985). In this solution, an initial overdensity causes a contin- 
uous sequence of spherical shells of cold matter, both pressure-free dark matter and baryonic 
fluid, centered on the overdensity, to slow their expansion, turn around and recollapse. The 
baryonic infall is halted by a strong accretion shock, while density caustics form in the 
collisionless dark matter, instead, due to shell-crossing. The postcollapse, virialized object 
we wish to model is then identified with the spherical region of shell-crossing dark matter 
and shock-bounded baryonic fluid in this infall solution (for ratio of specific heats 7 = 5/3) 
for which the mass is the same as that of our top-hat and the trajectory of its outermost 
spherical mass shell was identical to that of the outer boundary of our collapsing top-hat 
at all times until it encountered the shock. According to the detailed similarity solution for 
this infall problem, this spherical region of post-shock gas and shell-crossing dark matter, 
it turns out, is very close to hydrostatic and isothermal and has virtually the same radius 
as that of the minimum-energy solution for the matching TIS. This offers some support for 
our "minimum-energy" anzatz and explains the dynamical origin of the boundary pressure 
as that which results from thermalizing the kinetic energy of infall. 

With this "minimum-energy" anzatz, we found that top-hat perturbation collapse leads 
to a unique, nonsingular TIS, which yields a universal, self-similar density profile for the 
postcollapse equilibrium of virialized cosmological haloes. Our solution has a unique length 
scale and amplitude set by the top-hat mass and collapse epoch, with a density at every 
point which is proportional to the background density at that epoch. The density profiles 
for gas and dark matter are assumed to be the same (no bias). The final virialized halo has 
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a small but flat-density core. For the EdS case, this core radius r is about 1/30 of the total 
size r t of the halo (i.e. the truncation radius). [We note that this core radius r = r 0i King/3, 
where r ,Kin g is the "King radius" defined by Binney & Tremaine (1987, equ. [4-124b])]. 
The central density p is about 500 times the density p t at the surface and 18,000 times 
the background density at the collapse epoch. At intermediate radii, p drops faster than 
r~ 2 , as fast as r~ 2 - 5 . The one-dimensional velocity dispersion ay of the dark matter and 
the gas temperature T are given simply by a v = ksT ' jm = 47rGpo r o> where m is the mean 
gas mass per gas particle. This temperature is significantly different from that predicted 
by the standard uniform sphere (SUS) and singular isothermal sphere (SIS) approximations 
adopted previously in the literature, with T = 2.16 T sus = 0.72 T SIS . 

The derived mass profile of our TIS solution as a function of its virial temperature (or 
velocity dispersion) and its formation epoch reproduces to remarkably high accuracy (i.e. 
of order 1%) the cluster mass-radius-temperature relationships and average mass profile for 
CDM haloes which Evrard, Metzler, and Navarro (1996) derived empirically by fitting a large 
set of detailed numerical gas dynamical and N-body simulation results of cluster formation 
in the CDM model. The TIS and NFW halo mass profiles are also in very close agreement 
(fractional deviation of less than ~ 10%) at all radii outside of a few TIS core radii (i.e. 
about one King radius) for NFW concentration parameters 4 < c NF w < 7. In short, our TIS 
model is in good agreement with the average properties of simulated CDM haloes although 
it differs from these numerical results at very small radii. 

The purpose of this paper is to generalize the results of Paper I to the case of a low-density 
universe, either open and matter-dominated (f2 < 1, A = 0), or flat with cosmological 
constant (fi = 1 — A o)- 111 the former case, the internal structure of the TIS haloes will 
be the same as for the EdS results of Paper I, when radius and density are expressed in 
dimensionless terms in units of the core radius and central density. However, these latter 
dimensional quantities will be functions of the total mass and the formation epoch of the 
halo, functions which differ from the EdS results and depend upon the value of Q . These 
differences between the EdS solution for the TIS model in Paper I and that presented here 
for the open, matter-dominated universe (f2 < 1) can be understood as a reflection of the 
different solutions for the spherical top-hat problem in the two cases. For the case with a 
cosmological constant, however, the differences with respect to the EdS solution are much 
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more extensive. In this case, not only is it necessary to consider the differences which result 
from the different solution for the spherical top-hat problem, but it is also necessary to 
take proper account of nonzero A in the virial theorem, the conservation of energy, and the 
isothermal Lane-Emden equation, as well. In § 0, we shall summarize the spherical top-hat 
perturbation problem generalized to these low-density Friedmann universe cases. We will 
also briefly describe how the conservation of energy and the virial theorem are generalized 
in the presence of some uniform background component X of energy density, such as the 
cosmological constant or quintessence, and how these are combined with the top-hat solution 
to derive the properties of the postcollapse virialized object when the latter is assumed to 
be a uniform sphere [i.e. the "standard uniform sphere approximation" in the presence 
of the X-component] . We will then specialize these general results to the case in which the 
uniform component X is the cosmological constant. In § [| we will generalize the Lane-Emden 
equation for an isothermal sphere in hydrostatic equilibrium to take account of the presence 
of the X-component and specialize our discussion of its solutions to the cosmological constant 
case. In § |], we shall apply the virial theorem to these generalized isothermal Lane-Emden 
spheres, taking care to account properly for the important effect of finite boundary pressure. 
In § |5|, we will derive the minimum-energy TIS solutions for these low-density universe 
cases, including convenient analytical fitting formulae for practical application, in which all 
properties of the virialized post-collapse object are given as explicit functions of the mass and 
collapse redshift of the object for a given background cosmology. Our results and conclusions 
are summarized in § |||, where we compare the TIS solutions for low-density universes derived 
here with the SUS and SIS approximations for those well as with the TIS solution of 

Paper I for the EdS case. Finally, in view of the suggestive physical correspondence between 
our TIS solution in the EdS case and the shock- and caustic-bounded sphere in the self- 
similar, spherical infall solution of Bertschinger (1985), we have added an Appendix [X| in 
which we show how the latter solution (for an EdS universe) can also be applied at early 
times in the generalized low-density universe cases described here, by a simple rescaling of 
parameters. 
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Cosmological Haloes in a Low-Density Universe 7 

2 SPHERICAL TOP-HAT PERTURBATIONS IN A LOW-DENSITY 

FRIEDMANN UNIVERSE AND THE STANDARD UNIFORM SPHERE 
APPROXIMATION 

2.1 Before Collapse: the Exact Nonlinear Solution 

The spherical top-hat model, an uncompensated spherical perturbation of uniform overden- 
sity within a finite radius (Gunn & Gott 1972), affords considerable insight into the dynam- 
ics of the gravitational growth of cosmic structure, while still having an exact, analytical 
or semi-analytical solution (e.g. Peebles 1980). In what follows, we shall consider top-hat 
perturbations in cosmological models with two components: (1) a nonrelativistic component, 
which comprises all forms of matter, luminous or dark, baryonic or non-baryonic, that can 
cluster under the action of gravity, and (2) a uniform component, which does not clump 
at any scale of interest. As a special case we shall consider low-density matter-dominated 
universes with no uniform component, as well. 

The equation of state of the uniform component X relates its pressure px to its energy- 
density px according to 

Px = ~ 1) Pxc 2 = wp x c 2 . (1) 
Let us also define the quantity 

Px.efr = Px + 3p x /c 2 = (1 + 3w)p x = (n - 2)p X - (2) 

The mean rest mass energy density pb of the nonrelativistic component and the energy 
density px then each vary with time according to 

p b (t) oc a(t)- 3 , (3) 

Px (t) oc a(t)- n , (4) 

where a(t) is the Robertson- Walker scale factor, and n is a non-negative constant. [For 
details and references, see e.g. Martel & Shapiro (1998), and Wang & Steinhardt (1998)]. 
Particular values of n correspond to models with a nonzero cosmological constant (n = 0, 
w = —1), domain walls (n — 1, w — —2/3), string networks (n — 2, w — —1/3), vacuum 
stress or massive neutrinos (n — 3, w — 0), radiation background (n — 4, w — 4/3), and 
quintessence (0 < n < 3, — 1 < w < 0). 
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The time-evolution of the scale factor a(t) is described by the Friedmann equation. For 



the two- component models considered here, this equation takes the form 



H(ty = Hi 



;i - no - n X)0 ) (-) ~\n (-) ~ 3 + n x , (-) " 

Van/ Van/ Van/ 



(5) 



2 ^^2 3 

+n — 

^CLq/ \clq/ Van, 

where H(t) is the Hubble parameter, the density parameters of the nonrelativistic matter and 
X components are Q = pb/ Pait and fix = Px/Pcrit, respectively, where p cr i t = 3H 2 (t)/8nG, 
and all subscripts "0" refer to present values. The redshift z that corresponds to a scale 
factor a(t) is given by the usual relation 

1 + Z = J TTY ^ 
a{t) 

In what follows, it will be convenient to define the scale factor at present by 

'O \ V(3-n) 



a o = \ 1 , matter-dominated, n = 1 ; \J) 

i - n 

— — — , matter-dominated, n < 1 . 

n 

While the X-component does not clump under the force of gravity, its presence not only 
affects the the rate of expansion of the universe as indicated by equation (f|); it also modifies 
the gravitational forces on the matter component, as follows. The Poisson equation in the 
presence of this non-clumping component is 

V 2 $ tot = 47rG(p + p x ,efr), (8) 
and the corresponding gravitational potential is 

$tot = $ + y = $ + $x, (9) 

where $ is the gravitational potential due to the nonrelativistic matter component as it would 
be in the absence of the X-component, and $x is the correction due to the X-component. 
From Birkhoff's theorem, the collapse of a spherical top-hat density perturbation can be 
described by the Friedmann equation for a universe with a higher average density than 
that of the background universe outside the top-hat. Special cases of the top-hat model 
in universes with a nonclumping component were discussed in Lahav et al. (|1991|) (for the 



cosmological constant), and Wang & Steinhardt (1998) (for quintessence). Here we shall 
only give the basic equations and refer the reader to these papers for further details and 
references. 

The density inside the top-hat stays uniform during collapse. We can, therefore, describe 
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its evolution solely in terms of its overdensity 5 with respect to the background. The equation 
for the evolution of the top-hat overdensity in the presence of a nonclumping component 
with n < 3 is (Shapiro, Martel, and Iliev 2001) 

d 2 H6 + 1) _ 1 f dln(5+l) \ 2 + 9l dln(5+l) _ 3(5+1) = Q 

din a 2 3 V din a J gi din a g2 
where 

g x = 2na + 1 + (4 - n)a 3 "", <? 2 = 2(«a + 1 + a 3 ""), « = U ^ ^jff . (11) 

The boundary conditions are given by 

6(0) = 0, ^(0) = A, (12) 

where A indicates the initial amplitude of the top-hat fluctuation. For k—0, the solution 5(a) 
of equation (|i~0f ) is independent of f2o and Qx,o- In that case, we need solve equation (|i~0"D 
only once to obtain the family of solutions for different Qq. For open, matter-dominated 
models, equation fllOD is replaced by 

d 2 ln(5 + l) 1 fd ln(5 + l)V 2a + 1 d In (5 + 1) 35 

din a 2 3\ dlna J + 2 (a + 1) cilna 2(a + 1) ~ ^ ^ 

(Shapiro et al. 2001). 

Let us define the critical density contrast 5 CT i t as the solution 5l of the linearized version 
of equation ( |T0"D or equation (|T3|), extrapolated to the epoch at which the nonlinear solution 
5 predicts an infinite overdensity. This quantity 5 cr i t labels the time of collapse in scale-free 
units. The values obtained at turnaround will be denoted by subscript "m" for "maximum 
expansion," while the values at collapse time will be denoted by subscripts "coll" . Henceforth, 
we shall refer to z co \\ as the redshift which corresponds to the epoch of infinite collapse, at 
which 5 = oo, at time i co u. 

2.2 After Collapse: The Standard Uniform Sphere Approximation 

The properties of the final, virialized object which is postulated to result from top-hat 
collapse were derived in the "standard uniform sphere approximation" for the case with a 
cosmological constant by Lahav et al. (|1991|) and for the case of quintessence by Wang & 
Steinhardt ( |1998|) . The postcollapse state is fully described in this approximation by the 
radius and internal velocity dispersion of the final equilibrium sphere. 



© 0000 RAS, MNRAS 000, 000-000 



10 I.T. Iliev and P.R. Shapiro 

The SUS approximation assumes that the collapse of the top-hat to infinite density is 
interrupted by a rapid equilibration when 5l = o~ C iit, which results in another uniform sphere 
in virial equilibrium. The final radius r v - 1T of the virialized sphere is obtained by assuming 
that the energy of the top-hat is conserved during collapse and virialization and applying 
the virial theorem to the final state, assuming that the boundary pressure term in the virial 
theorem can be neglected. The conserved total energy of the sphere is E = K + W to t, where 
K = t/ t h + ^kin, U t h and T kin are the thermal and kinetic energy, respectively, 

Wtot = I [ pV$ tot • rdV = W + W x (14) 



2Jv' 

is the total potential energy in the presence of the X-component, the terms W and Wx are 



defined by the integral expression in equation ([14]) with $ and $x, respectively, and $tot 5 $ 
and $x are defined by equation ([J). 

Let us define the collapse factor which relates the size of the sphere at maximum ex- 
pansion r m to its postcollapse equilibrium size r vir as r\ = r vir /r m . The relative gravitational 
importance of the X-component and the matter component at maximum expansion are in- 
dicated by the dimensionless ratio 9 = — (px,eff/2p) m , where p is the matter density inside 
the top-hat. At the point of maximum expansion [which, unlike the EdS case, exists only 
for top-hat perturbations of amplitude high enough to enable collapse (see, for example, 
Martel 1994)], the sphere is cold and at rest, so K = 0, and its energy consists entirely of 
gravitational potential energy. For a uniform sphere of mass M and radius r m , equations (§) 
and (HD yield 



III 



After collapse, when the system settles down to a virial equilibrium, the total potential 
energy is 



vir 



1+077 



3 / Qcoll \ ' 

V a m J 



(16) 



5 r vir \ 2p J v . r 5 r v 
where a m is the scale factor at turnaround and a co u is the scale factor at collapse. According 
to the virial theorem, the kinetic and gravitational potential energies are related at this 
epoch by 

3( 7 -l)K + W -2W X = 0, (17) 
where 7 is the ratio of specific heats. We take 7 = 5/3. Together with the conservation of 
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energy, this implies that E = W/2 + 2Wx, and, therefore, the total energy in equation (|I5]) 
is related to the virial radius and densities after collapse according to 

^-ML^l . (18) 

Equating the total energy E in equations flT5|) and (|18D and using equation (f|) then 
yields an equation for the collapse factor 77, 

46n 3 (^) ~" - 2(1 + ^)77 + 1 = 0. (19) 

Hereafter we shall limit our discussion to the case n = (cosmological constant), with 
open, matter-dominated models as a limiting case for which px,cs = 0. Equation fllTf ) then 
becomes 

A6rf - 2(1 + % + 1 = 0, (20) 

where 9 = p\/p m and p\ is the constant vacuum energy density associated with the cosmo- 
logical constant [i.e. p\ = A/(8ttG)]. For any matter-dominated universe, 9 = and equation 
( p0|) reduces to the well-known relation 77 = 1/2. In general, equation ( pOf) has a closed form 
solution r/(^). However the expression is quite complicated and not very practical. Instead, 
simpler, approximate solutions have been proposed (Lahav et al. 1991, Kochanek 1995). We 
have derived our own, better approximation to the exact solution [with an error of order 

77 = 0.5 - 0.250 - O.1250 2 + O.1250 3 + O.218750 4 . (21) 

The upper limit for the 9 parameter is 1/2. If 9 > 1/2, then the right hand side of equa- 
tion (|8]) evaluated at the epoch of maximum expansion would be negative, and, therefore, 
the perturbation would not be bound. For 9 —> 1/2 from below, the top-hat perturbation is 
only marginally bound and would collapse arbitrarily far into the future. Plots of 9 vs. z co \\ 
for several different values of Ao = 1 — are shown in Fig. |l|. 

The virial temperature and velocity dispersion in this SUS approximation are derived as 
follows. The kinetic energy in the virialized state is 

K = E-W tot = |^9.[1 + 9(1 - 3 V % (22) 



m 



This kinetic energy is the energy of random internal motions only (T^ Q = 0). To obtain the 
equivalent temperature, we express the thermal energy as 
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Figure 1. Dimcnsionless parameter 9 = px/pm, the ratio of the cosmological constant energy density to matter rest-mass 
energy density inside the top-hat at maximum expansion, versus the collapse redshift of the top-hat, z co n, for different flat 
background universes, as labeled with the values of Aq: Ao = 0.6 (short-dashed line), Ao = 0.7 (solid curve), and Ao = 0.8 
(long-dashed line). 

U th = ^M , (23) 
2 m 

where m is the mean mass per gas particle. If run is the mass of a hydrogen atom, then 
m = fimH, where /i is the mean molecular weight. The virial temperature is, therefore, given 
by equations (|22|) and (p3| ) as 



T 



2 GM m 



[l + 0(l-3ry 2 )]. 



(24) 



5 k B r n 

This temperature is always higher than the corresponding temperature with no cosmological 
constant. The maximum departure from the case with A = is 30%, for 9 = 1/2, for which 
rj = 0.366. Henceforth, we shall refer to the virial temperature in equation (|24D as Tgus- For 
the case of a collisionless gas, we replace the virial temperature above by the virial velocity 
dispersion, 

(v 2 ) k B T 



4 



(25) 
3 m 

For a given 9, the mean overdensity of the virialized object in the SUS approximation 



with respect to the critical density of the universe at z co \\ is 



P 



Pcrit(^coll) 



077 3 (0) 



PcriH^collj 
Pcrit.O 



-1 



(26) 
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where a is defined by equation (|7|), according to which Qo^l = Ao if ^o + Ao = 1, and where 



equations 

Pcrit(^coll) 
Pcrit,0 



and (y 
h 



yield 



(1 - fio - A )(l + ^coll) 2 + fi (l + ^collf + A . 



(27) 



Equations (|20|) , (26) and (pT|) form a set of simultaneous equations which may be solved for 
the dependence of rj and 9 on z co \\, once A c is known by integrating the nonlinear top-hat 



perturbation differential equation (jlOf) [or ("p~3|)"| - In general, this solution for the dependence 
of r) and 9 on 2 co n for a given A c involves the roots of a cubic equation which can be 
cumbersome to express analytically. However, for 9 « 1 (z co ii >> 1), a good approximate 
solution is given by 
1 

2' 



rj = 2q 
and 



(28) 



rf 
where 



^7 



M^coll) 

h 



1 -2 



(29) 



(30) 



The dependence of the quantity A c on z co \\ can be expressed in terms of approximate ana- 
lytical fitting formulae (with errors ~ 1%) according to 

A c = 18vr 2 + c x x - c 2 x 2 , (31) 

where x = Q(z co n) — 1, and C\ = 82 (60) and c 2 = 39 (32) for the flat (open) cases, fi + Ao = 1 
(Q Q < 1, A = 0), respectively ( [Bryan fc Norman 1998| ), where 

^o(l + ^) 3 



n(z) 



(32) 



(1 - fio - A )(l + zf + O (l + zf + Ao' 
Henceforth, we shall refer to the quantities rj and A c for the SUS approximation as r/sus an d 

A c ,sus- 



3 ISOTHERMAL SPHERES IN THE PRESENCE OF THE 
X-COMPONENT 

For matter-dominated universes, both open and flat, the final virialized object decouples 
from the expanding cosmological background from which it condensed. Hence, when we 
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describe it as an isothermal sphere in hydrostatic equilibrium, we do so in the usual non- 
cosmological way (e.g. Binney and Tremaine 1987). Similarly, in the presence of an X- 
component, the final virialized object also decouples from the expanding cosmological back- 
ground from which it condensed, but it continues to be affected by the X-component be- 
cause of the modification that component causes to the gravity force, as discussed earlier. 
We can still describe it as an isothermal sphere in hydrostatic equilibrium in the usual non- 
cosmological way, but we must account properly for this modification of the gravity force. 
The Poisson equation in equation (|]) for the gravitational potential in the case of spherical 
symmetry is 

* =-47rG(p + p XiCfr ). (33) 



}_d_ / 2 ri$tot 
r 2 dr \ dr . 

The equation of hydrostatic equilibrium, Vp = pg, where g = — V$tot, combines with 
equation (|33|) in spherical symmetry to become 

k B Tdp GM(r) 4irGrpx, c s fQA , 

-j- = -pV$ tot = ~P 2 P 5 ' 34 

m dr r z 3 

where M(r) is the mass inside radius r. Multiplying equation ( |34] ) by r 2 m/ pksT and taking 
the derivative with respect to r, we obtain 

Let us consider the case of collisionless particles, too. The equilibrium velocity distribu- 
tion of the particles is a Maxwellian distribution given by 

ft \ Po /$tot-wV2\ , Q(n 

/w= (2^^ exp ( oj y (36) 

where p is the central density if we take $tot( r = 0) = 0, and oy is the one-dimensional 
velocity dispersion. After integrating over velocity, we obtain 

p = f f(v)dv = poe*™/°v, (37) 



which we substitute into equation ([33|), to obtain 
d ( 2 d{\np) \ 4vr 2 

Tr [r—)=-^G( P + Px , e{i )r. (38) 
By calculating the mean square velocity we obtain 

(v 2 ) = 3a 2 v . (39) 
The equivalent temperature can be calculated from 
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(v 2 ) 3k B T 



2 2 m 
to obtain for oy: 



(40) 



4 = — • (41) 

m 

A comparison of equation (^) with equation (^Sj) using equation ([II]) shows they are iden- 
tical. Hence, the structure of a self-gravitating isothermal fluid sphere in hydrostatic equi- 
librium for 7 = 5/3 is identical with that for a spherically-symmetric system of collisionless 
particles in virial equilibrium with three-dimensional particle orbits which are isotropic and 
have a spatially uniform velocity dispersion. 

To make equation (|38D nondimensional, we introduce new variables 

P = — , C = — > PX.eff = , (42) 

Po r Po 

where p is the central density, and we choose 

r = ^/(4vrGpo) 1/2 . (43) 
In terms of these variables, equation fl38|) becomes 



We must solve equation (f44|) with the following boundary conditions: 

P(0) = 1, ^(0) = 0. (45) 



An important difference between equation ( |4*4j ) and the standard isothermal Lane-Emden 
equation should be noted here. The standard form of the equation for the case without the 
X-component does not provide a natural cutoff of the density profile, and, thus, both the 
profile and the mass are infinite. In the presence of an X-component, however, the modified 
Lane-Emden equation shows that the density profile is finite as long as n < 2, since there 
is a matter density below which the matter cannot be gravitationally bound. According to 
equation flilD , this cutoff occurs at p(r) = —px, c s- For the special case of a cosmological 
constant, this cutoff occurs where p = 2p\. The mass of any such halo is, therefore, also 
finite and depends on the value of px,eff- 

In Figure 0, we specialize to the cosmological constant case and plot the dimensionless 
density profile p/po, the logarithmic slope of this density profile, and the dimensionless 
circular velocity profile v c /o~y, normalized to the velocity dispersion ay, all as functions of 
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12 3 4 
log 10 (r/r D ) 

Figure 2. Modified Isothermal Lane-Emden Spheres in the Presence of a Cosmological Constant, (top panel) Density profile (in 
units of the central density po) versus radius (in units of the core radius ro) when the X-component is a cosmological constant, 
for p\ = 0.1 (label 1), p\ = 0.001 (label 2), p\ = 10 — 5 (label 3), p\ = 10 — 7 (label 4), and p\ = (no cosmological constant; 
label 5). Each plot extends out to its natural cutoff at p = 2p\, except the matter-dominated one (curve 5), which is infinite and 
is plotted only out to r/ro = 10 4 . (middle panel) The logarithmic slopes of these density profiles versus dimensionless radius, 
(bottom panel) The circular velocity profiles v c (normalized by the velocity dispersion of the halo cry) versus dimensionless 
radius. Note the steep fall of the circular velocity close to the cutoff radius due to the fact that outer shells are only marginally 
bound in the presence of cosmological constant, unlike the case A = 0. 

r/r for several representative values of p\. In the presence of a cosmological constant, this 
circular velocity is: 



d$ GM(r) 



1 



air r \ p{r) J ^ ^ 

When px,cff = 0, there is a well-known, analytical, power-law solution to the isothermal 

Lane-Emden equation, if we relax the inner boundary conditions in equation (^) so as to 

permit a singularity at the origin. This singular isothermal sphere (SIS) solution is given by 



p{r) 



a 



v 



2ttGV 2 ' 
and 
2 1 GM(r) 



a 



in- 



constant. 



(47) 



(48) 



2 r 

For the modified isothermal Lane-Emden equation in the presence of the X-component, 
however, no such analytical solution is possible. Nevertheless, since p(r) << p\ over most 
of the profile, except for the outer parts, the SIS solution is still a reasonable approximate 
solution of equation ((44]) out to p ~ p\. As for the case px,eff = 0, this approximate SIS 
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solution when px.eff 7^ also does not satisfy the inner boundary condition in equations fl45*D . 
We will compare our results for the nonsingular case with this approximate SIS solution in 
§ |6.1| below, after we complete our derivation of the nonsingular TIS solution. 

As mentioned above, unlike the standard result when px, e s — 0, isothermal spheres in 
the presence of the X-component always have a finite radius and mass. Hence, unlike the 
case with no X-component, we do not need to truncate these spheres in order to apply 
them to describe realistic finite structures. However, as will be shown below, our final TIS 
model actually imposes a smaller truncation radius r t than that which corresponds to the 
marginally bound mass shell in the isothermal sphere solution. The total mass Mq of the 
isothermal sphere is then 

M = M(r t ) = P A7rp(r)r 2 dr = 4irp r%M(( t ), (49) 
Jo 

where ( t — r t/ r o and M(( t ) is the dimensionless mass: 

^^rt medC (50) 

If such a truncation radius r t does exist, this leads to the necessity of an external pressure 
to keep the system in equilibrium, which requires the inclusion of a surface pressure term in 
the virial theorem. This correction and its consequences are discussed in the next section. 



4 THE VIRIAL THEOREM FOR TRUNCATED ISOTHERMAL SPHERES 
IN A LOW DENSITY UNIVERSE 

Let us consider the general isothermal sphere density profile, p(r), obtained in the previous 
section. From the ideal gas law, the pressure inside the sphere as a function of the radius is 

p(r) = p{r) = o-yp(r), (51) 

and at the outer edge 

Pt = p(r t ) = o-yp(r t ). (52) 

The mean pressure inside the sphere is 

JpdV _ 3tfp(Q( 2 d( _3M(Ct) _ a{ (53) 

where ( t = r t/ r o- 



© 0000 RAS, MNRAS 000, 000-000 



18 I.T. Iliev and P.R. Shapiro 

The virial theorem for a static sphere in the presence of a surface pressure and of the 
X-component now becomes 

3( 7 - l)K + W + S p - 2W X = 0, (54) 

where K is once again just Uth, and S p is the surface pressure term. The thermal energy for 
a gas with a ratio of specific heats 7 = 5/3 is given by 

U th = 7 pdV = — = -aptV, 55 

7 — 17 7 — 1 2 

where V is the total volume. The surface term is equal to 

S p = - J pr ■ d S = -Wp t . (56) 

Hence, according to the virial theorem, the unmodified gravitational potential energy term 
W is 

a - 1 

W = -2 U th + 2W X , (57) 

a 

and the total energy is 
2 — a 

E = Uth + SW x . (58) 

a 

In order for the halo to be bound, the condition E < must be met, which requires 

a > 1 lw x - ( 59 ) 

The virial temperature of this generalized isothermal sphere is 

2a GM m 2 
r= 5(a-2) fe B r m [1 + ° {1 ~ * )] ' (60) 

which we shall henceforth refer to as 

For 9 = 0, the temperature in equation (|60|) reduces to T TIS derived in Paper I. In the 
presence of the X-component, this temperature is always slightly higher than the corre- 
sponding temperature for px = 0, however, but never by more than 5%, with the maximum 
reached for 9 = 0.5. Since a /(a — 2) > 1 for any a, the temperature Ttis is always higher 



than Tsusi the standard value for a uniform sphere shown in equation (|2~^). 

Just as we did in Paper I, we identify the virial radius r vir with the size of the truncated 
isothermal sphere (i.e. r vir = r t ). For comparisons with the results of dynamical calculations 

* Henceforth, as in Paper I, the notation "TIS" shall refer to a solution of the isothermal Lane-Emden equation (modified here 
to account for the presence of the X-component), with the nonsingular boundary conditions of equation (^jj) at the origin. 
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of the formation of such an equilibrium object, we should interpret r vir (and r t ) as the radius 
inside which hydrostatic equilibrium holds (e.g. Cole & Lacey 1996). Hereafter, we specialize 
the discussion to the case in which the X-component is the cosmological constant. 

5 CHOOSING A UNIQUE PROFILE: THE MINIMUM-ENERGY 
SOLUTION IN A LOW-DENSITY UNIVERSE 

5.1 The Nonsingular, Truncated Isothermal Sphere Formed by Top-Hat 
Collapse 

The virialized object which results from the collapse of a given top-hat density perturbation 
must have the mass of the top-hat before it collapsed and virialized. We assume that the 
total energy E is also conserved, including the effect the cosmological constant has on the 
potential energy, if present. Fixing the mass M and the energy E of the TIS which results 
from the collapse of the top-hat, however, is not enough to specify which member of the 
infinite family of solutions of the modified isothermal Lane-Emden equation which all share 
this mass and energy is chosen by the collapsing top-hat. In order to make a unique choice, 
some additional information is required, such as the value of the boundary pressure p t . As 
in Paper I, we shall assume that the correct boundary pressure p t is that value which makes 
the TIS with this total energy E and mass M correspond to the unique minimum-energy 
solution. In order to determine the minimum-energy TIS which results from a given top-hat 
collapse, we must take account of the fact that this top-hat collapse depends upon the value 
of the parameter 9. This parameter measures the relative importance of the vacuum energy 
density associated with the cosmological constant and the rest-mass energy density of the 
matter inside the top-hat at maximum expansion, which is different for different epochs of 
maximum expansion even for a given value of the cosmological constant. As a result, the 
internal structure of the TIS solution which results from a given top-hat collapse even when 
expressed in dimensionless form with radius in units of r and density in units of po will 
depend, not only upon the background cosmology parameters (i.e. Qq = 1 — Ao and H ), 
but also upon z co \\ for the top-hat. This makes the calculation of our minimum-energy TIS 
solution in this case more complicated than for the EdS case in Paper I. 

As in Paper I, this conservation of the energy E of the top-hat before and after its collapse 
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and virialization assumes that there is no extra prfV-work that needs to be taken into account 
due to the presence of the external boundary pressure pt which would otherwise alter the final 
total energy E compared with its initial value before collapse. This is appropriate for the case 
at hand, since the collapse prior to the epoch of virialization is that of a cold, pressure-free 
gas which collapses supersonically. As in the well-known, self-similar, spherical infall solution 
of Bertschinger (1985), the original energy is converted from potential energy at maximum 
expansion to a mixture of infall kinetic energy and potential energy during infall, with a 
negligible share of the energy going into compressional heating. In that solution, the infall 
is halted by a strong shock. At this shock, the kinetic energy of infall is converted primarily 
into the thermal energy of the shock-compressed gas, with only a small portion remaining as 
kinetic energy of subsonic, postshock infall. While the self-similar infall solution is strictly 
correct only for the EdS universe, we expect similar behaviour in low density universes, as 
well. However, since no general self-similar infall solution exists in these cases, the analogous 
infall solution for low-density universes must be studied by numerical means. Of course, in 
the limit of high redshift, the self-similar solution of the EdS case can be applied even to 
the case of a low-density universe if properly rescaled (see Appendix |A|). 

By analogy with this infall solution, therefore, we identify the boundary pressure p t 
in the case of our TIS solution, not as a fixed external pressure which acts on the top- 
hat boundary throughout its collapse and virialization, but rather as something like the 
instantaneous post-shock pressure in the infall solution. As such, it has a physical origin in 
the conversion of the original energy of the collapsing top-hat, itself, from potential energy 
at maximum expansion into kinetic energy of infall during collapse and, finally, into thermal 
energy of the post-shock gas, always conserving the original energy E of the top-hat. 

5.2 Finding the Unique Minimum-Energy Solution 

5.2.1 Flat Universe with Cosmological Constant 

As shown above, the truncated isothermal sphere solutions form a one-parameter family, 
described by ( t = r t /r - the truncation radius in units of the core radius. Specifying ( t , the 
total mass, and the total energy completely determines the solution. Alternatively, we can 
specify the mass, the total energy and the applied external pressure p t . 
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The conservation of energy equates the energy after collapse and virialization to the 

energy of the top-hat at maximum expansion, E(t m ) = E y \ T . This yields 

3GMi L +£ A = ^ UaMo _ lfa , Gpj/l , (61) 
5 r m \ p m J a 

where 

I t = P pr 4 dr. (62) 
Jo 

In dimensionless variables, equation ( |6~1~1) becomes 

(1 + %tis = -^^ + ^I#, (63) 
where 

I t = J o Ct p( 4 d(, (64) 

and we have used the fact that 

Pa _ P^Pm _ 3g7j|ig M t 
Po Pm Po Ct 

For a given external pressure pt, the total energy of the modified isothermal Lane-Emden 
sphere is given by 

E = -^( 1+ ») = -^( 1+ 9)!5H 
5 r m 5 r Ct 

= ^(4 7 r 3 G)^M 3/ \ 1/4 (l + ^)^(C i ), (66) 

where 

E (Ct) = —i/^- (67) 

We have plotted the dependence of the dimensionless total energy E on ( t in Figure |3] for 
several representative values of the parameter 9. In order to indicate the dependence of the 
size of the sphere on ( t , we can nondimensionalize the radius r t , according to 
rjE_ = 3 
GMq 5 

This definition of also corresponds to the familiar dimensionless energy parameter used 
in discussions of the stability of isothermal spheres. 

For a given mass M and external pressure p t , the solution is specified uniquely only 
if we can uniquely identify a special value of ( t , or equivalently, of E. Apparently, for any 
truncated isothermal sphere of mass M which is confined by a given external pressure p t , 

© 0000 RAS, MNRAS 000, 000-000 



-A^-t^h— (1 + %tis. (68) 



22 I.T. Iliev and P.R. Shapiro 



-0.6 - 



O 



1 - 



1.2 - 



1.4 




Figure 3. The dimensionless energy E(Ct) versu s di mensionless truncation radius ft = rt/ro for several representative values 
of # (see also Table hh, chosen as described in §5.3 From top to bottom: 9 = (matter-dominated models), 9 = 0.0118, 



9 = 0.0249, 9 = 0.0604, 9 = 0.123, and 9 = 0.5. The symbols indicate the minimum-energy points for each curve. 

there is a unique value of Q which minimizes the total energy E. As in Paper I, we shall 
make the reasonable anzatz that this minimum-energy solution is the unique TIS solution 
preferred in nature as the outcome of the virialization of the sphere in the presence of a fixed 
external pressure. We offered evidence to support this anzatz in Paper I. It is possible that 
this is a general result for any TIS formed by relaxation in the presence of a fixed external 
pressure, but we are only concerned here with spheres that evolve from cosmological initial 
conditions. 

For each value of 9 in the allowed range < 9 < 0.5, the minimum value of E as 
a function of Q for a given p t is found by minimizing the dimensionless energy E(( t ) in 
equation (p7[), with r/ Tig found by solving equation flBBp. Unlike the EdS case discussed in 
Paper I, the minimum-energy solution in a low-density universe is not universal (except 
for the special case A = 0), but instead depends upon the time of collapse, as implicitly 
parameterized by the quantity 9. This makes our calculation of the solution more complicated 
than in Paper I, since equation ([44] ) depends upon p\ and must, therefore, be solved in 
conjunction with equation ( |65"D to obtain a self-consistent solution and, subsequently, to 
find the minimum-energy point. As can be seen in Figure 0, the minimum of the energy 
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corresponds in each curve to a unique value of the dimensionless truncation radius ( t and 
becomes more pronounced when the cosmological constant contribution is more important. 
The value of (t at which this minimum energy occurs increases with increasing p\, from 
( t = 29.4 for matter-dominated models (9 = 0) to Q — 31.7 for 9 = 0.5, the maximum 
allowed value. The resulting values for ( t , a(( t ), Vtis, R = Po/Pt, Px, and M t versus 9 are 
shown in Figure |j. The point 9 = corresponds to the solution for matter-dominated models 
without cosmological constant (either open or flat), where these dimensionless quantities are 
identical to those for the EdS case obtained in Paper I. 

Thus far, we have derived the unique, dimensionless minimum-energy TIS solution for the 
postcollapse virialized object which results from top-hat collapse at any epoch for different 
cosmologies, as parameterized by the value of 9 = p\/ p m , the ratio of the vacuum energy 
density to the matter rest-mass energy density of the top-hat at its epoch of maximum 
expansion. In order to complete the dimensionless description of the TIS model for the 
postcollapse equilibrium object created by a top-hat which collapses at a given redshift z co \\, it 
is necessary to relate 9 to z co ii, according to the solution of the nonlinear top-hat perturbation 
equation ( ]T0| ) [or (jTB])]. This dependence of 9 on z co \\ for the top-hat perturbation was 
described in § [2] above. 

The parameters we have derived for the unique minimum-energy TIS solution for several 
representative values of the parameter 9, chosen so as to span the allowed range, are summa- 
rized in Table [l]. These values of the parameter 9 were chosen to correspond to the EdS case 
and open-matter-dominated cases (9 = 0), and to the flat, A-dominated model currently 
favoured by astronomical observations, with Ao = 0.7, for different collapse epochs: z co \\ = 1 
(9 = 0.01183), z coll = 0.5 (9 = 0.0249), and z coil = (9 = 0.0604), as well as for the extreme 
cases of A = 0.9, z co \\ = {9 = 0.123), and of either A — > 1 or else 1 + z coa — > (i.e. in 
the future) {9 = 0.5). If we adopt Ao < 0.9 as a reasonably conservative estimate of the 
possible range allowed by astronomical observation, then the maximum value of 9 of interest 
corresponds to that for which z co \\ = 0, i.e. 9(X = 0.9, z co \\ = 0) = 0.123. 

The radius r t of the TIS solution is encountered by the outer boundary of the collapsing 
top-hat at a time somewhat earlier than the time t co \\ of infinite density in the top-hat 
solution. We shall refer to the time at which the radius of the collapsing top-hat equals 
r t as tcross, to distinguish it from t co n. At time t cr0 ss, the top-hat exact solution yields an 
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9 8 

Figure 4. Dimensionless parameters of the TIS model in a flat universe with cosmological constant vs. the parameter 8, where 
6 = pxl pm , the ratio of the vacuum energy density to the matter rest-mass energy density at maximum expansion of the parent 
top-hat. 

overdensity 5 = 5 cross ~ 100 for EdS universe, while the extrapolated linear solution at 
this time yields d~L,cross ~ 1-56. For a flat, low-density universe, however, both 5 croS s and 
^L.cross exceed their values for the EdS case. For example, if z co \\ = 0, then 5 cro ss = 197(416) 
and Across = 1-96(2.52) for Ao = 1 — ^o = 0.7(0.9), respectively. For Ao — > 1 or else 
in the future for any value of A , both <5 cross — > oo and Across - ► oo. If we assume that 
the TIS forms instantaneously at t — t CTOSS , then its mean density corresponds to a mean 
overdensity 5 = 5 cmss when compared with the background density at Across, but the same 
TIS corresponds to a mean overdensity 6 = S co n, which is larger than this 5 CTOSS by about 
30%, when compared with the background density at z co \\. This value of 5 co \\ differs somewhat 
from the conventional value of 8(t co i\) found for the postcollapse virialized sphere in the SUS 
approximation. The values of <5 cross , <5L iCross , and S co u for our chosen set of illustrative values 
of 9 are also given in Table |l|. 

As noted in Paper I, applications of the top-hat model involving the Press-Schechter 
approximation customarily identify the characteristic time of formation of objects of a given 
mass as the finite time at which the nonlinear perturbation solution predicts collapse to 
infinite density, at which time the linear solution yields Sl = S CTit . Since SL jCross < S CTit , 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological Haloes in a Low-Density Universe 25 

it may be appropriate to replace £ cri t in such applications by the value 5L )Crossi implying 
that the number of objects formed at any epoch with a mass greater than some value may 
be somewhat higher for the TIS solution than previously assumed in applications of the 
Press-Schechter approximation. 

5.2.2 Open, Matter- Dominated Universe: Low-Density Universe without X-component 



For a matter-dominated universe, 9 = 0, and the family of solutions obtained in § |5.2.1 
reduces to a single dimensionless solution, which is independent of Qq, identical to the 
solution we obtained in Paper I. The dimensional solution, however, depends upon the 
background cosmology. The details of this dependence will be discussed in § [5~3| . As for 
the case of the flat, low-density universe discussed above, both S CTOSS and Across f° r the 
low-density, matter-dominated case exceed their values for the EdS case. For example, if 
^coii = 0, then 5 cross = 222 and Across = 3.29 for A = 0, Q = 0.3. As Q(z) — > (i.e. in the 
future), we again find that both <5 croS s —* oo and Across — > oo. 

5.2.3 Stability 

The TIS model in an EdS universe represents a stable solution of the isothermal Lane- 
Emden equation (Paper I). This conclusion holds for any matter-dominated model, as well. 
In the presence of an X-component, however, the Lane-Emden equation is modified, leading 
to a case whose stability, as far as we know, has not been studied. Nevertheless, since the 
TIS solution in the presence of the X-component departs from the solution for a matter- 
dominated universe by only a small amount for the observationally-constrained values of Ao, 
we expect our solution for the low-density universe cases considered here to be stable, as 
well. However, detailed study of this problem is beyond the scope of this paper. 

5.3 Dependence of the TIS Model Parameters on Halo Mass and Collapse 
Redshift for Different Background Cosmologies 

In this section, we describe how the dimensional parameters for our dimensionless TIS so- 
lution are specified for a given mass Mo and collapse redshift z co \\. For a given 6, the di- 
mensionless minimum-energy TIS solution and the nonlinear top-hat perturbation solution 
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combine to yield the mean overdensity of the TIS solution with respect to the critical density 
of the universe at z co \i, 



A 



c,TIS 



P 



Pcrit 



PcritV^collJ 
Pcrit.O 



(69) 



(where f2o a o = -^o f° r ^o + ^o = 1)- Note that this overdensity is somewhat different from the 
corresponding overdensity A Cj sus derived from the standard uniform sphere approximation 
as described in § pT2] , since the collapse factors rj = r m /r t are different for the two cases. 
These two overdensities are simply related, however, as follows: 



A _ ^SUS a 

^TIS 



(70) 



The dimensional parameters for the TIS solution for a given total halo mass Mo and A c .tis 
can be expressed as follows: 



3M 



1/3 



v4vrA CjT isPcrit(2 ; coii) 

n Ct 
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Tr 
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-a 



v c (r) 



k B 
GM(r) 



,_Pa_ 

' p{r) 



1/2 



(71) 

(72) 

(73) 

(74) 
(75) 

(76) 



For practical application of the TIS solution, we have provided a convenient set of accu- 
rate analytical fitting formulae for the dependences of the dimensionless parameters of the 
solution (e.g. ( t , Vtis, M t ) on 9 in Appendix JTJ. Since fitting formulae are also available for 
the dependence of A Cj sus on z co \\ for different cosmological parameters (Ao,^o), in the form 
of equation (|3~1~D, it is useful to express the dependence of 9 on z co a for a given background 
cosmology by combining equations (E6f) and (|2T|) to write 



9 



A c ,sus%us 



h(z co \i) 



h 



(77) 
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Table 1. Summary of the minimum-energy TIS solution in a flat universe with a cosmological constant for several illustrative 
values of 9. 



Quantity 


e = o 


6 = 0.0118 


6 = 0.0249 


e = 0.0604 


e = 0.123 


9 = 0.5 




29.400 


29.534 


29.677 


30.04 


30.58 


31.69 




3.730 


3.720 


3.709 


3.684 


3.649 


3.585 


J st/s 


2.156 


2.163 


2.170 


2.188 


2.213 


2.262 


M(Ct) 


61.48 


61.76 


62.06 


62.80 


63.90 


66.13 


PA 





1.405 x 10" 5 


2.842 x 10~ 5 


6.166 x 10~ 5 


1.042 x 10~ 4 


1.710 x 10~ 4 


7ns = 7^ 


0.5544 


0.5490 


0.5432 


0.5277 


0.5016 


0.3800 


VSXSS 


0.5 


0.4970 


0.4937 


0.4845 


0.4676 


0.3660 


n.Tis 

r t,SUS 


1.109 


1.105 


1.100 


1.089 


1.073 


1.038 


, cross 


1.562^| 


1.630 


1.710 


1.955 


2.524 


CO 


<5 (Across ) 


102.6 a 


116.2 


133.5 


196.8 


416.1 


CO 


<5tis(*co1i) 


129.6 


148.5 


172.7 


264.4 


579.1 


CO 


£* 

PO 


1 O/ifi v in - 3 
l.y4D X 1U 


1 OQ/1 v in _ 3 


1 oon v in _ 3 
i.yzu x iu 


1 Q07 w i n— 3 
l.ooi X IU 


l.OOl X IU 


1. 1 oy x iu 


R, 


514 


517 


521 


530 


544 


575 


PO 
Pb.coll 


1.796 x 10 4 


2.064 x 10 4 


2.430 x 10 4 


3.784 x 10 4 


8.641 x 10 4 


CO 


Pt 
Pb.coll 


35.0 


39.92 


46.67 


71.39 


158.8 


CO 




-0.3326 


-0.3333 


-0.3340 


-0.3357 


-0.3380 


-0.3420 


if/|W| 


0.6832 


0.6846 


0.6861 


0.6897 


0.6948 


0.7066 



a These values are slightly different from the values quoted in Paper I, Table 1. The correct values are shown here. If we match the 
outer radius of this top-hat at its turnaround epoch to the radius of the mass shell in the self-similar infall solution of Bertschinger 
(1985) which contains the same mass at the same turnaround epoch, which we did in Paper I, that mass shell crosses the shock at 
5 L = 1.5572. 

In that case, for a given z co n, fl , X , and h, equation ([77]) combines with equation flHf ) which 
relates ^sus to 9 and the fitting formulae for A Ci sus as a function of z co \\ in equation ( [311) to 
yield a pair of simultaneous algebraic equations for 9. 

For most cases of current interest (i.e. Qq — 0.3, which corresponds to 9 < 0.06 for 
-^coii > 0), a further approximation to these more accurate results is possible which is even 
simpler to use, if we neglect the effect of the cosmological constant on the dimensionless 
solution of the Lane-Emden equation. In particular, as long as we take proper account of 
the evolution of a top-hat density perturbation in a low-density universe up to its moment of 
infinite collapse and still satisfy energy conservation and the virial theorem, we will capture 
most of the dependence of the TIS solution on the background cosmology. The result can 
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be expressed as follows: 
' M 



1/3 



n = 

r = 
T = 

4 = 
Po = 
where 

-FOcoll 



337.7 
187.2 
6.367 



1O 12 M 0/ 

Mp ' 
1O 12 M 0/ 

Mp ' 
10 12 M 



1/3 



1/3 



7.843 x 10- 



0.59 



[F(^coii)]-^- 2/3 kpc, 
[i^coiOr 1 /^ 273 kpc, 
[F(^coii)]"^" 2/3 kpc, 

2/3 



M 



1O 12 M 

2/3 



- F(z coll )h 2 / 3 K, 

3/ 



1.098 x 10 4 (yq^) F(^coii)^ 2/3 km 2 s- 2 , 
1.799 x 10 4 [F(^ coll )] 3 pw = 3.382 x 10 - 25 [F(^ coll )] 3 /i 2 g/cm 3 . 



M^coll) 

h~ 



A. 



c.TIS 



Zcollj ^o) 



A c ,tis(A = 0) 



Q n A 



c,SUS 



.1/3 



coll; 



(78) 

(79) 

(80) 

(81) 

(82) 
(83) 

(84) 



^coii) 18vr 2 

For the EdS case, F = (1 + z C oii)> while for an open, matter-dominated universe and a flat 
universe with a cosmological constant, F — > f2 ^ 3 (1 + ^coii) at early times [i.e. a: — > 0]. Here 
fi is the mean molecular weight, where \l = 0.59 (1.22) for an ionized (neutral) gas of H and 
He with [ife]/[iy] = 0.08 by number. The relative accuracy of these formulae in comparison 
with the exact numerical results is plotted for the case Ao = 1 — ^o = 0.7 in Figure [5], which 
shows that the errors are always small, even for z co \\ = 0. 



6 SUMMARY AND DISCUSSION 

6.1 Comparison of the Minimum-Energy TIS Solution with the Standard 
Uniform Sphere and Singular Isothermal Sphere Approximations 

The temperature derived here for the minimum-energy TIS solution for the postcollapse 
sphere in virial and hydrostatic equilibrium which follows top-hat collapse is a factor of 
approximately two larger than the value previously derived for the SUS approximation, i.e. 
by satisfying energy conservation and the virial theorem for a postcollapse sphere of uniform 
density for which the surface pressure term is neglected, as summarized in § |^. In particular, 
if we write Ttis = brTsus, then bx — a. /(at — 2) ~ 2.2. The size of the TIS sphere, given 
in terms of the radius r m of the top-hat at maximum expansion according to r t = //Tis r m ; 
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Figure 5. Fractional deviation of the approximate analytical fits in equations ( fFsj ) — ( jig ) from the exact numerical solution 
for background cosmology with Aq = 1 — Ho = 0.7, as labelled. 

is actually not far from the size of the SUS sphere, r vir = ^sus r m ; as shown in Table |l|. 
However, in terms of the TIS core radius r (where, recall, our r ,Tis = ^King/3 as defined for 
conventional isothermal Lane-Emden spheres), the truncation radius is large, r t /r > 29 A, 
implying that the TIS is very far from uniform density. Despite this relatively small core 
radius, the TIS solution is also quite different from that of a singular isothermal sphere[j].For 
the latter, the ratio of the average density to that at the surface is 3, while we find a value 
of 3.6 to 3.7 for the nonsingular TIS. The truncation radius in the singular limit is only 
r t = (5/12)r m , smaller than the radius of the TIS. Similarly, the correction factor for the 
virial temperature of this SIS relative to that of the SUS is bx = 3, as opposed to the 
value for the TIS model which is 0t ~ 2. The central density of the TIS is > 1.8 x 10 4 times 
larger than the mean density of the background at the collapse epoch for the parent top-hat. 
Finally, the TIS solution predicts that the top-hat will virialize somewhat earlier than the 
nominal collapse time of the top-hat, since the outermost mass elements encounter a shock 

t As discussed in § |j| the singular isothermal sphere is not an exact solution of the Lane-Emden equation in the presence of a 
cosmological constant, but for the purpose of this comparison, we ignore the effect of A ^ on this equation but take account 
of the more important effects of A ^ on the top-hat evolution, on the conservation of energy of the top-hat before and after 
its collapse and virialization and on the virial theorem. 
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Table 2. A comparison of three approximations for the postcollapsc equilibrium structure of top-hat density perturbations in 
a low density universe. 

Uniform Singular TIS 



SphereQ Isothermal Solution 
Sphere[] 6 = - 0.5| 



^/^SUS 


1 


0.833 


1.11-1.04 


T 

Tsus 


1 


3 


2.1b— 2.2b 


po 


1 


oo 


514-575 


pt 








<£> 


1 


3 


3.73-3.59 


Pt 








n 


- NA - 


oo 


29.4-31.7 










A c 
A c,SUS 


1 


G) 3=L728 


1.36-1.12 


K/\W\ 


0.5 


0.75 


0.683-0.707 



a A top-hat perturbation of a given mass collapses at a given redshift in a background universe with given values of Qo and 
Ao; all of these values are held fixed in this comparison of the three approximations. 

b These SIS numbers are an approximation which ignores the small modification of the Lane-Emden equation solution to take 
account of A ^ 0, but accounts for the more important effects of A ^ on top-hat evolution, energy conservation and the virial 
theorem. 

c These values show the full range allowed for top-hat perturbations which collapse at finite time including the distant future 
(i.e. A -> 1). 

and shell-crossings in the infall solution at finite radius. This implies that the standard value 
of 5 cr jt ~ 1.69 used in the Press-Schechter approximation for the halo mass function at a 
given time, based upon extrapolating the linear growth to the epoch at which the nonlinear 
top-hat solution predicts infinite density, should perhaps be replaced by (<5L,cross)Tis > 1-56. 
A comparison of the TIS model with the SUS and SIS approximations is summarized in 
Table 

6.2 Summary 

We have generalized the TIS model, derived in Paper I for an Einstein-de Sitter universe, 
to the case of a low-density universe, either matter-dominated (Qq < 1, Ao = 0) or flat 
with cosmo logical constant (Qq + Ao = 1). The formalism we have presented can also be 
used to obtain the corresponding TIS solutions for other background cosmology models 
with a component of the energy density which remains homogeneous on scales relevant to 
the formation of virialized halos, such as quintessence. The halo density profile we have thus 
derived has a universal, time-invariant shape in the matter-dominated cases, when expressed 
in units of the central density po with radius in units of the core radius r . However, in the 
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presence of a cosmological constant, whose importance increases with time, the isothermal 
Lane-Emden equation is modified, and, as a result, this dimensionless TIS density profile is 
no longer time-invariant, but, instead, depends on the epoch of collapse. This dependence 
is relatively weak, except for the outer parts of the halo at late times. For example, for 
A = 0.7 and z co \\ = 0, the dimensionless radius of the halo r t /r is 30.04, larger than 
the value of 29.40 for a matter-dominated universe by about 2% (see Table p]). At later 
times (i.e. in the future) in the A-dominated universe, however, the departure from the 
universal, time-invariant shape of the matter-dominated case will be more pronounced, with 
the dimensionless radius reaching 31.69, or ~ 8% higher. For all low-density universe cases, 
including those which are matter-dominated, far more significant differences from the TIS 
model in an EdS universe, however, are found in the over-all amplitude of the density profile 
as measured relative to the mean density of the background. For example, the ratios of 
the average density inside the virial radius and of the central density to the cosmic mean 
background density at the epoch of collapse, are larger then their values for an EdS universe 
by more than a factor of 2, for the cases Qq = 0.3, Ao = and 1 — Qq = Ao = 0.7, for 
halos which collapse today, while these ratios grow arbitrarily large in the future, due to the 
decrease of Q(t) with time in both cases. 

The dimensional parameters of the TIS halo solution also depend significantly on the 
background cosmology. For example, a TIS halo of lO lo /i~ 1 M which collapses at z = in 
an EdS universe has a radius r t = 40.35 h~ l kpc, central density po = 3.38 x 10~ 25 /i 2 gcm" 3 , 
and velocity dispersion ay = 22.6kms _1 , while a halo of the same mass which collapses at 
the same redshift, but in a flat universe with Ao = 0.7, is ~ 18% larger, with a 33% lower 
central density, and ~ 9% lower velocity dispersion. The same halo in an open universe with 
Q = 0.3 would be ~ 14% larger than in the EdS case, with 38% lower central density, and 
6% lower velocity dispersion. 

Our results demonstrate that the presence of a cosmological constant can influence the 
internal structure of virialized haloes, particularly in their outer regions. This is true even if 
we use a conservative estimate of the contribution of the cosmological constant to the mean 
energy density of the universe, for haloes collapsing at the present. In principle, this effect 
grows in importance for haloes which are destined to form in the future. 

The TIS model has many characteristics in common with the haloes with uniform- density 
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cores which are expected to form from self-interacting dark matter (SIDM) (Burkert 2000, 
Dave et al. 2000, Firmani et al. 2000). Since our model describes the final equilibrium state 
of a collapsing halo and is based on a minimum-energy principle, it may transcend the 
details of the particle interactions which lead to this equilibrium. If so, then it is reasonable 
to expect that our TIS solution will be a useful approximation for the equilibrium structure 
of SIDM haloes. A detailed comparison between the TIS model haloes derived here and the 
predictions of halo formation in the SIDM model would, therefore, be of interest. 

The application of the TIS model presented here to the problem of halo formation in the 
CDM model, including further comparisons with N-body simulation results, will be described 
elsewhere. In one such application, the TIS model is shown to provide a good theoretical 
explanation for the observed rotation curves of dark matter - dominated galaxies and the 
correlation which has been reported between the maximum velocity on a given rotation 
curve and the galactocentric radius at which it occurs (|Iliev fc Shapiro 2001|). 
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APPENDIX A: APPLICATION OF THE SELF-SIMILAR, SPHERICAL 
INFALL SOLUTION AT EARLY TIMES IN A LOW-DENSITY UNIVERSE 

The self-similar cosmological infall solution of Bertschinger (1985) can be generalized to the 
case of a low-density universe at early times by simple scalings of the variables. At early times, 
Q(z) ph 1, and both the background universe and density fluctuations evolve approximately 
as in an EdS universe. Therefore, the self-similar infall solution will be approximately valid. 
Consider the spherical mass shell which is just turning around at some redshift z. Let the 
radius of this shell be r ta . According to the top-hat solution applied to this spherical shell, 
the mass m ta enclosed by this shell is just 



where pb{z) is the mean background matter density at redshift z and (1 + 5 m ) is the av- 
erage overdensity inside this sphere at its epoch of maximum expansion. At early times, 
for which the behaviour approaches that of an EdS universe, (1 + 8 m ) = 9tt 2 /16 and 
Pb(z) = Pcrit(z) = 3H 2 (z)/ (8ttG). However, for a low-density universe at high redshift, 




(Al) 





(A2) 



so 




(A3) 
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Therefore, if we fix the mass inside r ta at redshift z, then we must take a different value for 
the radius r ta for the EdS and low-density universes, according to equations ( |A1|) and ([A 3D, 
which yields 

rta(^)/r t a,Eds(2 ; ) = [Pcrit(^0, ^)/Pcrit,Eds(2;)] _1/3 = ^O^- ( A4 ) 

Consider now the mass within a sphere whose radius r at this redshift z is a fixed fraction 
Xb [in the notation of Bertschinger (1985), where we have added the subscript "B" to this 
dimensionless radius coordinate to distinguish it here from the unrelated quantities A and 
X(t) used elsewhere in this paper to refer to the cosmo logical constant in units of the critical 
density, as measured at the present and at time t, respectively] of the radius r ta , 

m(X B , z) = m ta (z) T^p- , (A5) 
1 + o m 

where M(Xb) is a dimensionless function of Xb only, for which M(l) = 1 + 5 m . According to 
equation (|A"5|), if we fix the mass m ta inside r ta at redshift z when comparing the EdS and 
low-density universe solutions, we must also have the same mass m(XB,z) inside any Xb in 
the two solutions. In that case, the radius r for any Xb must correspondingly relate to the 
radius rEds which encloses the same mass at the same redshift, according to 

r(X B )=r EdS (X B )n 1/3 . (A6) 

The age of a low-density universe at a given redshift is related to the age of an EdS universe 
at the same redshift according to the solutions of the Friedmann equation for these two 
cases. Equation ( |A~!| ) is easily integrated analytically to yield the result for high redshift, 

t(Q ,z)=t EdS (z)Qo 1/2 - (A7) 

This determines the scaling of velocities in the infall solution, as follows. At any Xb, the 
velocity scales as 

v oc - (A8) 
Equations (|£§) - Q58|) yield 

v(X B ) =v EdS (X B )nl /6 . (A9) 
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APPENDIX B: ANALYTICAL FITTING FORMULAE FOR THE 
DEPENDENCE OF THE TIS HALO DIMENSIONLESS PARAMETERS ON 
AND THE DEPENDENCE OF THE DENSITY PROFILE ON THE 
DIMENSIONLESS RADIUS C 

Analytical fitting formulae which closely approximate the numerical results for the dimen- 
sionless TIS parameters will make applications of our TIS model much more convenient. We 
have derived such fitting formulae for two overlapping intervals of 0, with different fractional 
errors in these two intervals, as follows. 

In the currently most relevant cosmological range, < < 0.123, the fits are given by 



Ct = 29.4003 + 11.46520 - 13.84280 2 - 12.84530 3 , (Bl) 

a = 3.7296 - 0.8660690 + 1.927420 2 - 1.7O3260 3 , (B2) 

Vtis = 0.554384 - 0.455290 + O.212580 2 + O.O21280 3 , (B3) 

TZ = 513.842 + 282.0310 - 193.6170 2 - 748.2540 3 , (B4) 

p x = 0.001232630 - O.OO3853450 2 + O.OO5798870 3 , (B5) 

M t = 61-485 + 23.88870 - 32.98540 2 - 14.O2720 3 . (B6) 



The relative errors of these fits are < 0.01% (except for the fit to p A , for which the relative 
error as — > is higher, since the exact solution has a slightly different slope at = from 
that of the fitting formula, p A — > a9, where a ps 0.001217. However, for < 0.005 where the 
above fit for p\ is not perfect, the departure from the result for the EdS case is negligible 
anyway. ) 

The fits in the full allowed range < < 0.5 are given by 



Ct = 29.3893 + 12.04740 - 2O.7520 2 + 11.7320 3 , (B7) 

a = 3.72913 - 0.8430680 + 1.71O170 2 - 1.2O9630 3 , (B8) 

VT1S = 0.554437 - 0.458040 + O.239990 2 - O.O42790 3 , (B9) 

K = 513.366 + 306.8120 - 478.8710 2 + 222.7260 3 , (B10) 

p x = 10 _3 (1. 11560 - 2.6O930 2 + 2.12420 3 ), (Bll) 

M t = 61.4683 + 24.7820 - 43.85780 2 + 25.88690 3 . (B12) 



The relative errors of these fits are < 0.1% (again with the exception of the fit for p A ). 
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In Paper I we obtained a simple analytical fit to the TIS density profile, as follows: 

«« = ?T?-Ff?' (B13) 

for C < Ct, with 

(A, a 2 , B, b 2 ) T is = (21.38, 9.08, 19.81, 14.62). (B14) 

This fit is accurate to within 3% inside the virial radius for Ao = 0. In the presence of a 
cosmological constant, the dimensionless profile varies, as described above. Nevertheless, for 
Ao < 0.9 and haloes that collapse by z = 0, the fit still has accuracy of better than 3%. 
Overall the fit deteriorates a bit, however, by slightly overestimating p at all radii. 

The corresponding fit to the TIS halo circular velocity profile (i.e. the rotation curve), 
obtained simply by integrating equation ( |B13|) , is 

v e (r) = cry | A - B + — bB arctan (^r) - aA arctan (— ) } ' (B15) 

( |Llicv fc Shapiro 200 1| ). This fit to v c {r)jay is good to < 1% for A = 0, and for r < 2r t /3 
if Ao 7^ 0, while at larger radii, depending on Ao and z co u, the error could be larger but is 
not greater than 6% at r t , even for A = 0.9 and z co u = 0. Most of the departure of the 
exact circular velocity profile from this fit is due to the pA-dependent correction factor in v c , 
which was omitted here for simplicity. In Iliev & Shapiro (2001), we apply the TIS model to 
explain the observed rotation curves of dark matter-dominated galaxies and the statistical 
correlations amongst their rotation curve parameters for different mass haloes collapsing at 
different epochs in the CDM model for different background cosmologies. 
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